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ABSTRACT 

The stellar surface mass density profiles at the centers of typical ~ and lower-mass 
spheroids exhibit power law "cusps" with S oc R~ n , where 0.5 < 77 < 1 for radii ~ 1 — 100 pc. 
Observations and theory support models in which these cusps are formed by dissipative gas 
inflows and nuclear starbursts in gas-rich mergers. At these comparatively large radii, stellar 
relaxation is unlikely to account for, or strongly modify, the cuspy stellar profiles. We argue 
that the power-law surface density profiles observed are a natural consequence of the gravita- 
tional instabilities that dominate angular momentum transport in the gravitational potential of 
a central massive black hole. The dominant mode at these radii is an m = 1 lopsided/eccentric 
disk instability, in which stars torquing the gas can drive rapid inflow and accretion. Such a 
mode first generically appears at large radii and propagates inwards by exciting eccentricities 
at smaller and smaller radii, where (< R) -C Mbh- When the stellar surface density profile 
is comparatively shallow with 77 < 1/2, the modes cannot efficiently propagate to R = and 
so gas piles up and star formation steepens the profile. But if the profile is steeper than 77 = 1, 
the inwards propagation of eccentricity is strongly damped, suppressing inflow and bringing 
r] down again. Together these results produce an equilibrium slope of 1/2 < 77 < 1 in the 
potential of the central black hole. These physical arguments are supported by nonlinear nu- 
merical simulations of gas inflow in galactic nuclei. Together, these results naturally explain 
the observed stellar density profiles of "cusp" elliptical galaxies. 

Key words: galaxies: active — galaxies: evolution — quasars: general — galaxies: nuclei — 
galaxies: bulges — cosmology: theory 



1 INTRODUCTION 

Observations have established that typical < L» ellipticals and 
bulges exhibit steep central "cusps" in their surface luminosity den- 
sity and stellar mass density profiles - i.e. a continued rise in a 



power-law like fashion towards small radii dLauer et al.|1991| j l992 
Crane et al.||1993[ [Ferrarese et al.||1994| |Kormendy et al.|l994 



Lauer et al.| 1 995 ||Kormendy 1 1 999) , |Faber et al. | ( | 1 997| > showed that 



power-law nuclear profile ellipticals also tend to have higher de- 
grees of rotational support and diskyness. This, together with other 
observations jKormendy|1999||Quillen et al.|2000||Rest et al.|2001| 
|Lauer et al.|2007||Ferrarese et al.|2006||C6te et al.|2007| >, has sup- 
ported the idea that the cusp ellipticals are the direct product of 
gas-rich mergers and nuclear star formation during such mergers. 
Quantitatively, power-law cusps have 



lot E* oc/T 



(1) 



with 0.5 < 77 < 1 representing the typical observed slopes; the 
power-law profile extends from the smallest radii observed in 
nearby spheroids lpc) to anywhere from ~ 10 to ~ lOOpc 
(Ferrarese et al. |l994||Kormendy et al.|2009) nThe most massive 
spheroids deviate from this behavior and exhibit flattened nuclear 
profiles, or "cores." This is, however, widely believed to be due to 
"scouring" by a binary black hole in a gas-poor environment (see 
e.g. |Begelman et al.|19 80) and thus does not reflect the initial for- 
mation history of the central stars that we focus on here. 



* E-mail : phopkins ©astro . berkeley.edu 

1 Provided it is defined over the same dynamic range, this is non-parametric 
and the same logarithmic slopes are recovered when Sersic fits are used. 



|Barnes & Hernquist| ( |1991[ > and |Mihos "& Hernquist J1994| > 
showed in simulations that tidal torques in mergers can drive rapid 
gas inflows, providing the fuel to power intense nuclear starbursts 
and build up the central stellar surface densities (see also Kor- 
|mendy|[T99"9"l |Hopkins et aL|2009"l >. This is observed in local Ul- 
traluminous Infrared Galaxies (ULIRGs), whose nuclei constitute 
the most rapidly star-forming environments in the local Universe. 
Moreover, the observed central gas densities and star formation 
rates in ULIRGs will leave them with typical power-law like cusps 
when the starburst is complete (Hibbard & Yun 1999 ; Tacco ni et al.| 
|2002l >. 

Despite this progress, no theoretical explanation exists for why 
spheroid cusps should have a power law-like form in the range ob- 
served. The scales of the observed cusps are comparable to, or less 
than, the black hole (BH) radius of influence, and the potential is 
thus quasi-Keplerian. At these radii, stars almost certainly formed 
primarily dissipatively in a gas-rich disk, rather than via violent re- 
laxation of a pre-existing stellar disk. On the very smallest scales, 
two-body relaxation is expected to drive the system to a Bahcall & 
Wolf (1976) cusp; however outside of ~ lpc the relaxation time 
is 3> fHubbie- Moreover, the fact that the observed central cusps are 
disky and elliptical, often with significant rotational support (see 
references above) suggests that two-body relaxation has not had a 
large effect. Instead, an understanding of the observed stellar cusps 
appears to require the combined effects of angular momentum re- 
distribution and star formation in disky, gas+stellar+BH systems. 

Recently, we have shown that the formation of lopsided, ec- 
centric disks within the BH radius of influence is a ubiquitous fea- 
ture in hydrodynamic simulations of massive gas inflows in galax- 
ies; such lopsided disks lead to efficient angular momentum trans- 
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Figure 1. Top: Nuclear stellar mass profiles (cusps) produced in hydrody- 
namic simulations of gas inflow driven by gravitational instabilities in dis- 
sipative starbursts; a lopsided/eccentric disk mode dominates at < lOOpc. 
Simulations of different resolution/spatial scale are each shown only over 
the modeled dynamic range; the highest resolution simulations extend to 
~ 100 pc while more typical simulations begin at ~ 10 pc. The simulations, 
taken from Hopkins & Quataert ( 20 lOai, have a wide range of initial mass 
profiles, gas fractions, and stellar feedback models. However, in all cases 
the simulated mass profiles are consistent with cusps having power law- 
like slopes in the range 0.5 < r) < 1 (S oc R~ v ; dashed line shows these 
two cases, which roughly bracket the simulations). Bottom: Distribution of 
cusp slopes. The slopes are measured in each simulation at lOpc (fitting a 
power law to each from 3 — 30 pc or the maximum dynamic range allowed 
within this interval), and from the observations in Lauer et al. l'2007j at ap- 
proximately the same radius in both volume-limited and magnitude-limited 
sub-samples (showing the "cusp," as opposed to "core" component of the 
observed bimodal distribution). 



fer from the ga s to the stars, powering BH a ccretion rates of up to 



lOMQyr (Hopkins & Quataert 



20 10a I. Moreover, the stellar 



relics of these disks are reasonably similar to the nuclear disks ob- 
served on < lOpc scales around nearby supermassive BHs (Hop- 
|kins & Qua taert 2010 b}, particul arly the well-studied case at the 
center of M3 1 |Lauer et al. 1993). There are also many candidate 
nuclear disks observed in other syst ems ^L atiCT^alJ1996] 



Houghton et al. 2006] |TTiatte et al.||2000; Debattista et al. 



2005 



Afanasiev & Sil'chenko 2002; Seth et al. 2010; Ledo et al. 



2006 
|2010) . 



Figure[T](top panel) shows the stellar surface density profiles 
at the end of the "nuclear-scale" and "ultra-high" resolution sim- 
ulations of Hopkins & Quataert (2010a), which extend inwards 
from > 10— 100 pc with ~ 0.1 pc resolution; we show results in the 
quasi steady state phase of all simulations with significant inflows, 



> 0.3Mq yr~ into < 1 pc, sustained for > 10 yr. These SPH sim- 
ulations include gas, stars, star formation, and a black hole as an 
additional collisionless particle; the simulations are idealized prob- 
lems focused on studying the nonlinear evolution of gravitationally 
unstable systems in the potential of a massive black hole. Hopkins 
|& Quataert]j201l) show that the central dynamics and inflows are 
dominated by the nuclear m = 1 modes. In Figure [T] the absolute 
stellar mass densities depend on the initial conditions (e.g., total gas 
mass), but the slopes are more robust; the simulations shown span a 
wide range in initial gas fractions, prescriptions for star formation 
and gas physics, initial stellar and gas mass profiles, and bulge-to- 
disk ratios (see Tables 1-3 in Hopkins & Quataert 2010a ), but con- 
verge to similar slopes. Comparing with the observed power law 
slopes of ellipticals (bottom panel), the agreement is reasonable. In 
this Letter, we provide a physical explanation for these results. 

2 PROPAGATION OF INSTABILITIES 

Physically, the lopsided or eccentric disk mode (azimuthal 
wavenumber m = 1 or amplitude oc cos <j>) is unique in any nearly 
Keplerian potential (Tremaine 2001). Gravitational torques from 
other modes are suppressed by the gravity of the BH. However, 
the resonant response between the epicyclic and orbital frequen- 
cies allows for global, low frequency m=\ modes that can exert 
strong torques on the gas by inducing orbit crossing and shocks 
(e.g., |Chang et al1|2"007l |Hopkins & Quataert|[2fjTT) . Because of 
the importance of the m = 1 modes for redistributing gas inside the 
potential of the BH, we now focus on the physics of these m = 1 
modes, in particular their propagation to smaller radii. 

2.1 The WKB Limit 

Consider an initially axisymmetric, thin, planar disk (surface den- 
sity E) with a BH of mass Mbh at the coordinate center; we use 
cylindrical coordinates throughout (R, (j>, z). The initial potential in 
the disk plane can be written $o = &o(R), and other properties are 
defined in standard terms: 



V: =R 



9$ _ GM enc (<fl) 
dR~ R 



■ Vc/R 



1 L dyn 

dR 



or 2 + 



(2) 
(3) 
(4) 



where V c is the circular velocity, f2 the angular velocity, and k 
the epicyclic frequency. We use c s to denote the sound speed in 
a gaseous disk and a- the vertical dispersion in a stellar disk. 

We consider a linear perturbation E — > Eo(/?) + Et(/?, 4>) 
(where E is the total gas+stellar disk surface density) in a frame 
rotating with the perturbation pattern speed Q p , and decompose the 
perturbation into linearly independent modes: 



E„, = E«(i?) exp{i(m(j) — cdt)} 



E a (i?) = \a(R)\ Eo(tf) exp<^ i / kdR 



(5) 



(6) 



where m is the azimuthal wavenumber, \a\ — \a(R)\ the effective 
mode amplitude, k the radial wavenumber, and the complex w the 
mode frequency. With these definitions, the mode pattern speed is 
Q p = Re(uj)/m, and the mode growth rate 7 = Im(oj). 

We adopt a power-law disk as a convenient reference model: 



E oc R 71 = E 



(7) 
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It is straightforward to show then that 



Q 2 



GM BH 27tqGE ( R 



Ro 



Ro 



-(i+i) 



(8) 



where a = (r[l - f ] T[^l]) /(F[^] r[f ]) for < 77 < 2. 

We first consider modes in the WKB limit of tight-winding 
(i.e. local modes), where \kR\ 2> m. We caution that this limit 
does not, in fact, hold for many of the global modes of most 
interest, but it is nevertheless instructive. We follow |Tremainel 
( |2001} 's derivation for a slow mode (Q p <C fi) in which the non- 
Keplerian part of the potential is small, i.e. $ = $bh + &d where 
Qd/^BH ~ Mj/Mbh <C 1. Expanding the equations of motion in 
terms C($rf/$BH) gives the WKB dispersion relation (to leading 
order in |/t/?| - ') of quasi-Keplerian slow modes, 



UJ = ZU + TT GT,J \k 



Q, — ct k z CI 



for a gas disk, or 

LU = VJ + TT G Erf I k\ T 

^vj + 7vGE d \k\Q' [ exp(-/3|fctf|) 
for a stellar disk, where we define 

2 -k 2 1 (2 d d 2 



(9) 



(10) 



20 



20, \r dr dr 2 



In the dispersion relation for a stellar disk (eq.| 10|>, T is the standard 
reduction factor (Binney & Tremaine 1987 1, and the latter equality 
in equation [TO] is a convenient approximation for softened gravity, 
with ft m cr z /Vc ~ h/R (the stellar disk scale height). 

The m = 1 slow modes are stable in the limit Mj <C Mbh 
(Tremaine 2001). Because of this physical constraint, the m — 
1 modes first appear at large radii - the radius = j^ cl it where 
Mrf/MsH ~ 1, i.e. where the potential is transitioning to Keplerian 
(Hopk ins & Quataert|2010a||Hopkins|2010} . The pattern speed O p 
of the unstable mode is ~ f2(i? C rit)- But if the mode can propagate 
inwards at constant fl p , it will eventually be a slow mode, relative 
to the local O at smaller radii. 

How does this propagation occur? The wave packets propa- 
gate with approximate group velocity ~ V c (.f?crit), so the timescale 
for the mode to travel is just the dynamical time at R a \ t - However, 
if these modes are forming in realistic "initial" disks, and if they 
are the dominant source of angular momentum transport, then the 
initial disk surface density profile cannot already be steep. Dimen- 
sionally, vj ~ {\/RO)d$ e /dR ~ GE/OR (exactly true if the disk 
has a locally power-law profile; vj = —a (2 — rj) irCE/OR). Since 
O diverges oc r~ 3/l2 at small radii, if the surface density profile is 
sufficiently shallow (and the dispersion is finite) then the right-hand 
side of Equation[9]becomes arbitrarily small as r — ¥ 0, and finite u> 
cannot be supported - the wave will refract back at some minimum 
radius R m \ n - From Equation|9] this g-barrier occurs when 



\0„\ >zu + 



n 2 (GS rf ) 2 

4 c 2 n 



in gas or 



\flp\ >w + 



ttGE 
e/3RO 



'■ vj-\- - 



ttGE 



(11) 



(12) 



in stars. Since zu ~ GT,/OR oc R l ^ 2 ~ TI at small radii, for systems 
with finite c s or constant /3 and a shallow E oc R~ v with 77 < 1 /2, 
the initial waves cannot reach R — 0. Note that Ostriker et al.H1992) 
show that the same restriction applies for modes in a pure fluid disk 
with a hard outer edge. 



If the slope is too shallow, but the m = 1 modes are present, 
they will drive gaseous inflows that will "pile up" near the refrac- 
tion radius. This will steepen the mass profile and increase the self- 
gravity at this radius, eventually allowing further mode propagation 
(in both gas and stars). Once 77 > 1 /2, then the RHS of Equation|9] 
no longer vanishes as r — ► 0, and the modes can propagate through 
to R = 0. Physically, the propagation can be understood as the ec- 
centric mode at larger radii exciting strong eccentric perturbations 
at smaller radii. 

Consider two nearly-adjacent annuli at radii R' and Ri : the ma- 
terial at i?i is part of the m = 1 mode, the material at < R' remains 
unperturbed. In the WKB limit the mode behavior at larger radii is 
swamped by the nearest asymmetric term - i.e. just inside Ri, the 
perturbing potential is w = 2nGT,i(R[) \k\~ l (since there 

is no local corrugation to cancel this out). In the global limit the 
result is similar: at small radii inside an eccentric ring at radius Ri 
having mass M r j ng and m = 1 amplitude \a\, the magnitude of the 
local perturbed potential is just w \a\ GMn„ g /Ri ~ 7rG£i(i?i)l?i, 
i.e., the same as the WKB result with \k\ ~ R^ 1 ■ To estimate the 
velocity induced at smaller radii by this perturbed potential, we 
note that for a cold gas or stellar disk, the local pattern speed of 
the m = 1 mode is just u — w ~ 7rGEo \k\ fi -1 . Together, this 
leads to the result that the response (in both gas and stars) at 
smaller radii ~ R' to the eccentric disk at larger radii is given 
by: \v r /V c \ ~ (Ei/E )|/ti?| _1 ~ \a\. In other words, for a non- 
negligible mode amplitude \a\ ~ E1/E0 and a global mode with 
\kR\ ~ 1, large eccentricities and hence large coherent m = 1 mode 
amplitudes, can be induced. The induced modes at these somewhat 
smaller radii can, in turn, excite large coherent eccentricities in the 
material at yet smaller radii, and so on, allowing the perturbation to 
grow even at arbitrarily small R. 

The above derivation also implies, however, that there is 
a regime in which the inwards propagation of eccentricity will 
be inefficient. For a global mode at R\, the perturbed potential 
is ~ tyGT,,(Ri)R u so the response |e| oc Ei(/?i)i?i/E(i?)i? = 
\a(Ri)\'£(Ri)Rif£(R)R. For a sufficiently flat mass profile 
E(Ri)Ri >T,(R)R,i.e.,r)< 1 for E(R) (xR~ v , the induced pertur- 
bation is large down to arbitrarily small R. But if the mass profile is 
too steep, E(/?i)i?i <~E(R)R, or 77 > 1, then although the mode can 
formally be supported even at small radii, the induced amplitude 
will decline as one moves to R — > 0. Crudely, we expect the prop- 
agation efficiency defined as log(|a(/J)|/|a(/?i)|) to decline oc —77, 
for 77 > 1 . This is explicitly demonstrated for a very large sample of 
models by Zakamska & Tremaine ( 2004 1 in the context of eccen- 
tricity propagation from an external perturber to the inner planets 
in planetary systems (the cutoff at 77 < 1/2 is not evident in this 
approach, however, because of the discrete nature of the problem). 
If 77 > 1 is established by some means, the low efficiency of ec- 
centricity propagation to small radii implies that gas will not inflow 
as efficiently at small radii. This will flatten the gas density profile 
and star formation will do the same for the stellar density profile, 
providing a mechanism for the system to self-adjust to have 77 < 1. 

2.2 Linear Global Modes 

The WKB results above are not exact, especially when the modes 
of interest are global and the disk mass is significant relative to 
Mbh (both of which are typically the case!). To show that our 
conclusions are robust, we also demonstrate the same points re- 
garding the propagation of modes using exact linear solutions for 
particular global normal modes. Our methodology is described in 
detail in Hopkins (2010), which we briefly summarize here. We 
define Ro for the power-law disk model (eq. |7J so that Mj(< 



© 0000 RAS, MNRAS 000, 000-000 



4 Hopkins and Quataert 



Rq) = Msn/la (2 - 77)] - i.e. E = Mbh/ '(2tt olR 2 ). Unlike in the 
WBK analysis we do not expand the equations to linear order in 
Md/MBH but keep the full linear perturbation equations. We con- 
sider a stellar-dominated (collisionless) disk since in our numerical 
results, the dominant torques in the gas are due to the stellar modes. 
The resulting equations of motion for linear perturbations are 



+ ■ 



o=-(n-w)s, 
s 

E 

7a 

E 



[2f2 {vt. + va - v&) - (fi- w)] $ fl 
(n - w) (1 + - va) + n (2 + - 



20 



1 a 

(13) 



where = dr' r' P(r, r') E a (V) follows from Poisson's equa- 
tion^ = d\nX/dlnR, and A = k 2 - (Q - uj) 2 . 

It is straightforward to solve eq. [13] for the eigenfunctions 
(normal modes) of the system. For convenience and realism we 
modify the disk mass profiles with a steep outer power-law cut- 
off (E oc R~ v (1 + [R/a] 2 )- (3 - v)/2 ) so that they have finite total 
mass Md — Mbh; Hopkins (2010) shows that the exact choice of 
Md and/or the cutoff radius has no affect on any of our conclusions. 

Figure[2]shows some of the resulting normal modes for a stel- 
lar disk. The growth rates 7 and pattern speeds Q,, are indicated 
on each panel, in units of Q(Ro). For any choice of disk parame- 
ters, there is a large variety of normal modes; here, we focus on 
the most rapidly growing "global" modes in each model. The re- 
sults are qualitatively similar for all global modes (local modes, 
potentially supported at all radii but localized in radius, are not of 
interest here). We take /3 = 0.1 for the softening, but our conclu- 
sions are essentially identical for a wide range of /3; in |Hopkins| 
( |20 10| > we show that this extends to /3 > 0.3, i.e. nearly-spherical 
configurations. This is because the manner in which E(i?) enters 
the equations means that the important dimensional parameter is re- 
ally M enc (< R) at a given radius, so puffier systems and even multi- 
ple overlapping disks making a quasi-spherical configuration give a 
qualitatively identical result. The important parameter we focus on 
here is the power-law index of the disk mass profile 77, for which we 
show various choices in Figure[2] 77 = 0.05, 0.50, 0.80, 1.10, 1.85. 

For each value of 77, Figure[2| shows the absolute value and real 
component of the surface density perturbation a(R) = E n /E and 
the induced eccentricity R a /R, where R a is the magnitude of the 
radial perturbation from the linear equations of motion. The modes 
are normalized so that MAX(|a(U)|) == 1. Where the eccentricities 
are significant, there can be orbit crossings and shocks in the gas. 
This dissipation helps drive rapid inflow; see [Hopkins & Quataert] 
( |2011) for a detailed discussion of this physics. 

The key result in Figure [2] is how the structure of the global 
modes changes with 77; this confirms our intuition derived from the 
WKB approximation. When the disk surface density profile is shal- 
low (77 < 1/2), the modes cannot propagate inwards efficiently - 



2 The kernel P is denned by P(r.r') = —it Gb l /i( r < I r >) / r > + 
7r Gr/r 1 2 and includes the direct and indirect components of the potential, 
respectively. The Laplace coefficient b t / 2 15 given by 

cos6d9 



2 r 

T Jo 



(14) 



(l-2xcos6+x 2 +j3 2 y/ 2 

where (3 represents the gravitational softening, as it appeared in the WKB 
approximation. For /3 > 0, the disk potential and Q are slightly modified 
accordingly, but this is a minor effect. 
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Figure 2. Global linear normal modes for gravitational perturbations to a 
nuclear stellar disk around a BH; the total disk-to-BH mass ratio = 1. Each 
row highlights one mode for a mass profile slope S <x R~ v and scale height 
ft ~ h/R = 0. 1 ; the pattern speed Q p and growth rate 7 of the mode are la- 
beled on the right in units of f2(/?o). Radii are in units of Rq, the BH radius 
of influence, ~ 10 — lOOpc. Left: Mode amplitude: black is absolute value, 
blue is Re(S /S). Right: Real and absolute value of the induced eccen- 
tricity. Shallow surface density profiles (77 < 1/2) cannot support modes at 
small R and so are spatially localized. For larger r\, the modes propagate to 
R — > 0. For r/ > 1 , however, the inwards propagation of eccentricity is less 
efficient (the disk is "stift'er" against external perturbations), which would 
suppress shocks and gas inflow at small R. The characteristic profile shapes 
can therefore self-adjust to have 1/2 < rj < 1. 



they are confined to a moderate range of radii. At 77 = 1/2, the 
modes are suddenly able propagate to arbitrarily small R. Going to 
somewhat larger 77 = 0.7 — 0.8, the structure of the modes is quite 
similar. For 77 > 1 , the induced eccentricity is strongly suppressed at 
small radii (see §|2Tj even though the mode formally has non-zero 
amplitude to arbitrarily small R. As a consequence, the induced gas 
inflow would also be strongly suppressed. 

3 DISCUSSION AND CONCLUSIONS 

The general interpretation of the stellar density profiles of < L„ el- 
lipticals is that violent relaxation produces the outer "wings" of the 
mass profile, while a nuclear starburst similar to that observed in 
nearby ULIRGs and many high-redshift galaxies produces a dense 
stellar relic that dominates the mass profile inside the central -~kpc 
(e.g., |Faberetal.|1997| l. The mass profile at smaller radii is thus set 
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by the physics of angular momentum redistribution and star forma- 
tion. This inner dissipative region is no longer sensitive to large- 
scale torques such as are produced in a merger. Instead, the inflow 
is likely due to secondary gravitational instabilities that develop as 
the disky component of gas and stars becomes self-gravitating on 
small scales (Shlosman et al. 1989|. 

In previous work (Hopkins & Quataert 2010aJ we have 
demonstrated that inside a radius from one to several times the BH 
radius of influence, the character of the instabilities that dominate 
angular momentum redistribution changes: the non-axisymmetry is 
dominated by an eccentric/lopsided disk or one-armed spiral mode 
- a "slow" m = 1 mode, unique to quasi-Keplerian potentials. 

If the stellar and gaseous density profiles are relatively shal- 
low (e.g., as might be the case absent an earlier epoch of gas in- 
flow) the m = 1 modes cannot be supported down to R — > 0, but 
reflect off of a g-boundary at a finite radius. Between this inner 
boundary and co-rotation, however, the modes will drive accretion, 
steepening the gas density profile; star formation will steepen the 
corresponding stellar density profile. As the mass profile steepens, 
the m = 1 modes can propagate deeper in the potential, until a crit- 
ical slope is reached, at which point the modes can propagate to, 
and drive inflow to, R = 0; for a power-law disk with E oc R~ v , the 
critical slope is r\ = 1/2. 

If gas is driven to small radii very efficiently, star formation 
will likely ensure that both the gas and stellar mass profiles further 
steepen. The surface density profile can, however, eventually be- 
come sufficiently steep that the inwards propagation of eccentricity 
is inefficient: the outer asymmetric perturbation is weak compared 
to the local disk self-gravity at small radii. For a power-law disk 
with E oc R~ v , this occurs at 77 = 1. The resulting pile-up of mass 
at larger radii, together with continued star formation, will flatten 
the gas and stellar mass profiles. 

The net result is a plausible equilibrium: in the presence of 
significant gas inflow from larger radii, star formation and the prop- 
agation of m = 1 modes will self-adjust so that the surface density 
profile satisfies 1/2 < r\ < 1 inside the potential of the central black 
hole. Remarkably, this range of slopes is comparable to what is ob- 
served in the centers of "cusp" ellipticals (Fig.[T](. 

On the smallest scales near the BH, dynamical relaxation 
plays an important role in setting the stellar density profile (see 
|Bahcall &~W olf 1976). Such effects are, however, unlikely to be 
important outside of ~ 1 pc because the A'-body relaxation time be- 
comes long compared to the Hubble time. The presence of massive 
perturbers can significantly accelerate stellar relaxation, but it is un- 
clear how disturbed the nuclei of these (now) gas-poor galaxies are. 
Moreover, the observational evidence for disky structures in "cusp" 
galaxies suggests weak/incomplete relaxation, which would have 
to be very extreme to operate out to 10 — lOOpc in any case. Even 
a merger of binary black holes tends to leave these structures if it is 
gas-rich (Hopkins & Quataert 2010a), because most of the nuclear 
gas inflows tend to follow the BH binary coalescence, regenerat- 
ing the disk, and multiple overlapping generations of disks can be 
formed (which will given the same mass profile, since each individ- 
ually must satisfy the same E (R) or Mdisk (R) scaling). And our con- 
clusions are robust to a number of variations in e.g. the disk thick- 
ness and detailed structure, In the potential of the BH, relaxation 
can be enhanced by the very same resonance between orbital and 
epicyclic motion that is so critical for the presence of m = 1 modes 
(Kocsis & Tremaine 2010). Scalar resonant relaxation, which mod- 
ifies the eccentricity axes of stellar orbits, is, however, inefficient at 
the radii of interest. Vector resonant relaxation is likely to be impor- 
tant, but this only changes the inclination angles of the disky orbits; 



this may wash out some of the observable diskyness or introduce 
warps into the nuclear kinematics, but it will not significantly affect 
the mass profiles. 

The combined effects of star formation and gas inflow driven 
by m = 1 modes provide a plausible explanation for the stellar mass 
profiles of "cusp" ellipticals at < 10 — 100 pc. This explanation 
links such cuspy profiles to the physics of angular momentum trans- 
port and BH growth. It also makes observational predictions: (1) 
the characteristic radii of those slopes should be correlated with the 
radii of influence of the BH (or radii enclosing comparable mass); 
(2) non-negligible radial anisotropy should be difficult to remove 
from the stellar orbits; (3) similar profile shapes should be observ- 
able in late-stage merger remnants, provided sufficient new stars 
have formed in the central regions. 
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